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Abstract 

In this paper we will present a quantum algorithm which works 
very efficiently in case of multiple matches within the search space 
and in the case of few matches, the algorithm performs classically. 
This allows us to propose a hybrid quantum search engine that inte- 
grates Grover's algorithm and the proposed algorithm here to have 
general performance better that any pure classical or quantum search 
algorithm. 



1 Introduction 

Quantum computers jHl El 112] are probabilistic devices, which promise to do 
some types of computation more powerfully than classical computers jHl E] ■ 
Many quantum algorithms have been presented recently, for example, Shor 
[T^ presented a quantum algorithm for factorising a composite integer into 
its prime factors in polynomial time. Grover ^U] presented an algorithm for 
searching unstructured list of N items with quadratic speed-up over algo- 
rithms run on classical computers. 
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Grover's algorithm inspired many researchers, including this work, to 
try to analyze and/or generalize his algorithm jH d 13 El- Grover's 
algorithm is proved to be optimal for a single match within the search space, 
although the number of iterations required by the algorithm increases; i.e. 
the problem becomes harder, as the number of matches exceeds half the 
number of items in the search space ^2] which is undesired behaviour for a 
search algorithm since the problem is expected to be easier. In this paper 
we will present a fast quantum algorithm, which can find a match among 
multiple matches within the search space after few iterations faster than any 
classical or quantum algorithm although for small number of matches the 
algorithm behaves classically. 

This leads us to proposing a hybrid search engine that includes Grover's 
algorithm and the algorithm proposed here. We also discuss the conditions 
that allow both algorithms to be integrate into a single hybrid. 

The plan of the paper is as follows: Section 2 gives a short introduction to 
quantum computers. Section 3 introduces the search problem and Grover's 
algorithm performance. Sections 4 to 6 introduce the proposed algorithm 
with analysis on its performance and behaviour. And we will end up with a 
conclusion in section 7. 

2 Quantum Computers 
2.1 Quantum Bits 

In classical computers, a bit is considered as the basic unit for information 
processing; a bit can carry one value at a time (either or 1). In quantum 
computers, the analogue of the bit is the quantum bit {qubit ^3]), which 
has two possible states encoded as |0) and |1); where the notation | ) is 
called Dirac Notation and is considered as the standard notation of states 
in quantum mechanics j7j. For quantum computing purposes, the states |0) 
and |1) can be considered as the classical bit values and 1 respectively. 
An important difference between a classical bit and a qubit is that the qubit 
can exist in a linear superposition of both states (|0) and |1)) at the same 
time and this gives the hope that quantum computers can do computation 
simultaneously {Quantum Parallelism). If we consider a quantum register 
with n qubits all in superposition, then any operation applied on this register 
will be applied on the 2" states representing the superposition simultaneously. 
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2.2 Quantum Measurements 



To read information from a quantum register (quantum system), we must 
apply a measurement on that register which will result in a projection of 
the states of the system to a subspace of the state space compatible with the 
values being measured. For example, consider a two-qubit system |</)) defined 
as follows: 



10) = a|00)+/5|01)+7|10)+5|ll), (1) 

where a, /3, 7, and 6 are complex numbers called the amplitudes of the system 
and satisfy \af + \l3f + \'-ff + \6f = 1. The probability that the first qubit 
of 10) to be |0) is equal to + If for some reasons we need to have 

the value |0) in the first qubit after any measurement, we must try some how 
to increase its probability before applying the measurement. Note that, the 
new state after applying measurement must be re-normalized so the total 
probability is still 1. 



2.3 Quantum Gates 

In general, quantum algorithms can be understood as follows: Apply a series 
of transformations (gates) then apply the measurement to get the desired 
result with high probability. According to the laws of quantum mechanics 
and to keep the reversibility condition required in quantum computation, 
the evolution of the state of the quantum system of size n by time t is 
described by a matrix U of dimension 2" x 2"^ 



1^') = U\ij), (2) 

where U satisfies the unitary condition: WU = I, where W denotes the 
complex conjugate transpose of U and / is the identity matrix. For example, 
the X gate {NOT gate) is a single qubit gate (single input/output) similar 
in its effect to the classical NOT gate. It inverts the state |0) to the state 
|1) and visa versa. It's 2x2 unitary matrix takes this form. 



X 



1 

1 



(3) 



and its circuit takes the form shown in Fig.(P). Notice that, from now on we 
assume that a horizontal line used in a quantum circuit represents a qubit 
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and the flow of the circuit logic is from left to right. For circuits with multiple 
qubits, qubits will be arranged according to the notation used in the flgure. 



{a\0)+P\l)) 



X 



iP\0)+a\l)) 



Figure 1: NOT gate quantum circuit. 



Another important example is the Hadamard gate {H gate) which has 
no classical equivalent; it produces a completely random output with equal 
probabilities of the output to be |0) or |1) on any measurements. It's 2x2 
unitary matrix takes this form, 



H 



1 
72 



(4) 



and its circuit takes the form shown in Fig.( 



x) 



Figure 2: Hadamard gate quantum circuit, where x is any Boolean variable. 

Controlled operations are considered as the heart of quantum computing 
[2], the Controlled-f/ gate is the general case for any controlled gate with one 
or more control qubit(s) as shown in Fig.(jHla). It works as follows: If any of 
the control qubits |cj)'s (1 < z < n — 1) is set to 0, then the quantum gate 
U will not be applied on target qubit |t); i.e. U is applied on \t) if and only 
if all |cj)'s are set to 1. The states of the qubits after applying the gate will 
be transformed according to the following rule: 

\ci) \ci) ;1 <i <n-l , . 

\t) ^ \tcu) = U'^'^-'-^ \t) ^''> 

where CiC2...c„_i in the exponent of U means the AND-'mg of the qubits 

Ci, C2, Cn-l- 



H 



^(io) + (-irii)) 



4 



\ci) - 

\C2) - 

I Cn-l)- 
1^) - 



u 



ci) 


\ci) 




\C2) 


Cn-l) 


Cn-l 


tcu) 


\t) 



|ci) 

|C2) 
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a.Controlled-U 



b.Controlled-NOT 



Figure 3: Controlled gates where the back circle • indicates the control 
qubits, and the symbol © in part (b.) indicates the target qubit. 

If U in the general ControUed-f/ gate is replaced with the X gate men- 
tioned above, the resulting gate is called a Controlled- A^OT gate (shown in 
Fig.(j21b)). It works as follows: It inverts the target qubit if and only if all 
the control qubits are set to 1. Thus the qubits of the system ci, C2, Cn-i, t 
will be transformed according to the following rule: 



\ci) — ^ |cj) ; 1 < i < n — 1 

1^) \tcN) = |t©CiC2...C„_i) 

where C1C2 . . . Cn-i is the AND-ing of the qubits Ci, C2, . 
classical XOR operation. 



(6) 

. , Cn-l and © is the 



3 Search Problem 

Consider a list L of N items; L = {0, 1, — 1}, and consider a function 
/ which maps the items in L to either or 1 according to some properties 
these items shall satisfy; i.e. f : L ^ {Ojl}- The problem is to find any 
i & L such that f{i) = 1 assuming that such i must exist in the list. It was 
shown classically that we need approximately N/2 tests to get a result with 
probability at least one-half. Let M denotes the number of matches within 
the search space such that 1 < M < N and for simplicity and without loss 
of generality we can assume that = 2". Grover's algorithm was shown to 

solve this problem in O {^^/WJM^ . In jlH], it was shown that the number 

of iterations will increase for M > N/2 which is un desired behaviour for a 
search algorithm. To overcome this problem it was proposed in ^3] that the 
search space can be doubled so the number of matches is always less than 
half the search space and then iterate the algorithm 7r/4-^2iV/M times so 
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Figure 4: Quantum circuit for the proposed algorithm. 



the algorithm still runs in O i^-^/NjMy But using this approach will double 

the cost of space/time requirement. In the following section we will present 
an algorithm that can find a solution for M > N/2 with probability at least 
92.6% after applying the algorithm once. 

4 The Algorithm 

4.1 Iterating the algorithm once 



For a list of size = 2", the steps of the algorithm can be understood 
as follows as shown in Fig.Q: 

1- Register Preparation. Prepare a quantum register oi n + 1 qubits all in 
state |0), where the extra qubit is used as a workspace for evaluating 
the oracle Uf-. 

\Wo) = |0)®" ® |0) . (7) 

2- Register Initialization. Apply Hadamard gate on each of the first n 
qubits in parallel, so they contain the 2" states, where i is the integer 
representation of items in the list: 
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Wr) = (i/®" ® I) \Wq) = ( ^ J] |z) ] ® |0) ; TV = 2". (8) 



i=0 



3- Applying Oracle. Apply the oracle Uf to map the items in the list to 
either or 1 simultaneously and store the result in the extra workspace 
qubit: 



m = Uf m = ^ 5^ (K) ® io © fm = ^ E (i^) ® 



j=0 1=0 



(9) 



4- Completing Superposition and Changing Sign. Apply Hadamard gate 
on the workspace qubit. This will extend the superposition for the n + 1 
qubits with the amplitudes of the desired states with negative sign as 
follows: 



N-l 



^ E (K) ® (|o) + |i))); P = 2N = 2«+^ 



VP 



i=0 



Let M be the number of matches, which makes the oracle [// evaluate to 
1 (solutions); such that 1 < M < A^; assume that Ej indicates a sum 
over all i which are desired matches (2M states), and indicates a 
sum over all i which are undesired items in the list. So, IPV3) can be 
re-written as follows: 



i=0 
N-l 

+^ E "(K)®(|o) + |i))) 

i=0 (11) 
N-l N-l \-^-^J 



;^E'(K)®|o))-;^E'(K)®|i)) 

i=0 i=0 
N-l N-l 



1=0 i=0 
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From Eqn. ljllj) : we can see that there are M states with amphtude 
(~Vvp) where f{i) = 1, and (P — M) states with amphtude (Vv^). 
Notice that, applying Hadamard gate on the extra qubit sphts the \i) 
states (solution states), to M states (K) ® with positive am- 
plitude (Vx/p) and M states ® |1))) with negative amplitude 

5- Inversion About the Mean. Apply the Diffusion Operator D similar to 
that used in Grover's algorithm ^Qj on the n + 1 qubits. The diagonal 
representation of the diffusion operator D can take this form: 

D = (2 |0) (0| - J) = 2 (iPl - I. (12) 

where, = -^J2k=o\^) equally weighted superposition of 

states. The effect of applying D [TH] on a general state "^^Iq ctk \ k) 
produces ^^Jq^ [— ctfc + 2 (a)] \k), where, (a) = J2k=o is the mean 
of the amplitudes of all states in the superposition; i.e. the amplitudes 
Uk will be transformed according to the following relation: 

Ofc ^ [-afc + 2 (a)] . (13) 

In our case, there are M states with amplitude (~Vvp) and P — M 
states with amplitude (Vv^), so the mean (a) is as follows: 

So, applying D on the system \W^) shown in Eqn.()ll|) can be under- 
stood as follows: 

a- The M negative sign amplitudes (solutions): will be transformed 
from (-1/^) to a , where a is calculated as follows: Substitute 
ak = ^ and {a) shown (Eqn.((T3))) in Eqn.lfT^ we get: 
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b- The (P — M) positive sign amplitudes will be transformed from 
(Vvp) to b , where b is calculated as follows: Substitute = 
and (a) shown fEqn. (jl4p ) in Eqn. (jl3p we get: 



~ vpy^ Pi- 



rn 



We can see that a> h after applying D. The new system \Wi) can be 
written as follows: 



N-l N-1 

D \W,) = \W,) = hY.' (K) ® |0)) + a E ' (K) ® 

N-1 N-1 I 

+&E"(K)®|o)) + 6E"(l^)®|i))- 

such that, 



Ma^ + {P- My = 1. (18) 

Notice that, if no matches exist within the superposition (i.e. M = 0), 
then all the amplitudes will have positive sign and then applying the 
diffusion operator D will not change the amplitudes of the states as 
follows: Substituting ak = and (a) = ^ (^P (^^^^ Eqn. lfT^ we 
get: 

^ ,i(pmi.4.„., (19, 



VP P\ \VpJJ Vp 

6- Measurement. Measure the first n qubits, we get the desired solution 
with probability given below: 



i- Probability Pg to find a match out of the M possible matches; 
taking into account that a solution \i) occurs twice as: (|i) ® |0)) 
with amplitude b and ® |1)) with amplitude a as shown in 
Eqn. ()17|l . can be calculated as follows: 
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= ^(l0-16(#)+8(ff) (20) 
= 5(f)-8(f)' + 4(fr. 

ii- Probability P„s to find undesired result out of the states can be 
calculated as follows: 

Pns = {P-2M)b\ (21) 

Notice that, using Eqn. |T8|) 

Ps + Pns = M{a^ + b^) + (P - 2M)b^ 

= Ma" + {P- M)b^ = 1. ^^^^ 

4.1.1 Performance after Iterating the Algorithm Once 



n, where N = 2"^ 


Max. prob. 


Min. prob. 


Avg. prob. 


2 


1.0 


0.8125 


0.875 


3 


1.0 


0.507812 


0.937500 


4 


1.0 


0.282227 


0.968750 


5 


1.0 


0.148560 


0.984375 


6 


1.0 


0.076187 


0.992187 



Table 1: Algorithm performance with different size search space. 



Considering Eqn. (fTH|) . Eqn. (fTH|) . Eqn. (pn|) and Eqn. (f?T|) . we can see that the 
probability to find a solution varies according to the number of matches M 
in the superposition. 

From Tabled we can see that the maximum probability is always 1.0, 
and the minimum probability (worst case) decreases as the size of the list 
increases, which is expected for small M because the number of states will 
increase and the probability shall distribute over more states while the aver- 
age probability increases as the size of the list increases. It implies that the 
average performance of the algorithm to find a solution increases as the size 
of the list increases. 
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To verify these results, taking into account that the oracle Uf is taken as 
a black box, we can define the average probability of success of the algorithm; 
average{Ps) , as follows: 



N 

average{Ps) = jn Yl ^CmPs 

M=l 

N 

— 2^ M\(N-My. -^^^ I," -t- o ; ^23) 

N 

where ^Cm = M\{N-My. number of possible cases for M matches. 

We can see that as the size of the list increases {N oo), average{Ps) shown 
in Eqn.(|2Sl) tends to 1. 

Classically, we can try to find a random guess of the item, which repre- 
sents the solution (one trial guess), we may succeed to find a solution with 
probability pj^'™') = M/N. The average probability can be calculated as 
follows: 



M=l 

N.M 



2N Yl 

M=l 
1 

2' 



^ (24) 

AI\{N-M)\N 



It means that we have an average probability one-half to find or not to 
find a solution by a single random guess even with the increase in the number 
of matches. 

Similarly, Grover's algorithm has an average probabilty one-half after 
arbitrary number of iterations as we will see. It was shown in j3] that the 
probability of success of Grover's algorithm after q iterations is given by: 



= sin2((2g + 1)^), where, < ^ < - and sm'^{9) = — . (25) 

The average probability of success of Grover's algorithm after arbitrary 
number of iterations can be calculated as follows (Appendix A) : 
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averageiPf) = ^Y1 ''CMsm\{2t + 1)^) = -. (26) 

M=l 

Comparing the performance of the proposed algorithm, first iteration 
of Grover's algorithm and the classical guess technique, Fig.(j3J shows the 
probability of success of the three algorithms just mentioned as a function of 
the ratio (M/N). 




M/N 



Figure 5: A plot of the probability of success of the proposed algorithm Pg, 
first iteration of Grover's algorithm Pf^'^' and the classical guess pj^'^''^***™') 
ClS Sb function of the ratio {M/N). 

We can see from Fig.® that the probability of success of the proposed 
quantum algorithm is always above that of the classical guess technique. 
Grover's algorithm solves the case where M = N/4 with certainty and the 
proposed algorithm solves the case where M = N/2 with certainty. The 
probability of success of Grover's algorithm will start to go below one-half 
for M > N/2 while the probability of success of the proposed algorithm 
will stay more reliable with propabilty at least 92.6%. For M < N/8, the 
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Figure 6: Quantum circuit for the iterative version of the proposed algorithm. 



probabihty of success of the proposed algorithm will start to go below one- 
half where performance of Grover's algorithm will be much better as we will 
verify in the next section. 



4.2 Iterating the algorithm 



If we consider iterating the above algorithm: For a list of size A^(= 2"), 
prepare n qubits and append extra q qubits for applying q iterations of the 
algorithm. The iterating version of the algorithm works as follows (as shown 
in Fig.®): 

1- Initialize the whole n + q qubits system to the state |0). 

2- Apply Hadamard gate on each of the first n qubits in parallel. 

3- Iterate the following, for iteration k: 

a. Apply the oracle U f taking the first n qubits as control qubits and 
the k^^ qubit workspace as the target qubit exclusively. 

b. Apply Hadamard gate on the k*^ qubit workspace. 

c. Apply diffusion operator on the whole n + k qubit system inclu- 
sively. 

4- Apply measurement on the first n qubits. 
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To understand how the iterative version of the algorithm affects the sys- 
tem, we will trace the state of the system during the first few iterations. 

Consider the system after the first iteration shown in Eqn. ()17|) . second 
iteration will modify the system as follows (to clear ambiguity, a and h used 
in the above section will be denoted clS QjQ and h^Q^ respectively, where the 
superscript index denotes the iteration and the subscript index is used to 
distinguish amphtudes) : 

1- Append second qubit workspace to the system: 



N-l 



N-l 



wn = bl^^ E ' (K) ® |o)) ® |o) + E ' (K) ® ® |o) 

' i=0 1=0 

+bS^ (K) ® |o)) ® |o) + bi;^ Yl" (1^) ® ® |o). 

i=0 i=0 



(27) 



2- Apply Uf as shown in step 3-a: 



N-l 



N-l 



= 6« E ' (K) ® |0)) ® |1) + a« E ' (K) ® ® |1) 

' j=0 j=0 

+b^^ E" (K) ® |o)) ® |o) + b^^^ (1^) ® ® |o). 

i=0 i=0 



(28) 



3- Apply Hadamard gate on second qubit workspace (J®"+i (g) H): 



W. 



(2) 



(1) iV-1 (1) iV-1 

= ^E'(K)®|o))®|o)-^E'(N)®|o))®|i) 

i=0 i=0 
(1) Af-1 (1) N-\ 

+^ E ' (N) ® ® |o) - ^ E ' (K) ® ® |i) 

j=0 4=0 

.(1) ^1 .(1) ^-1 



+^ E " (N) ® |o)) ® |o) + ^ E " (K) ® |o)) ® |i) 

i=0 'i=0 
(1) Af-1 (1) N-\ 

+^ E " (N) ® ® |o) + ^ E " (K) ® ® 



i=0 



i=0 



(29) 
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4- Apply diffusion operator as shown in step 3-c: 



+4'^ e'' (K) ® ® io) + af ) ' (I^) ® ® |1) 

i=0 j=0 
N-1 

+6?) E " (N) ® |o)) ® |o) + fei') E " (N) ® |o)) ® |i) 

+b'i^ e'" (N) ® ® io) + fef e' " (K) ® ® 

1=0 i=0 

(30) 

where the mean of the amphtudes to be used in the diffusion operator 
is calculated as follows: 



(«2) = ^ ((2"+2 - 4M) ^) 



(31) 



V2 A^^ 



And the new amplitudes Oq^'', af\h^Q^ and are calculated as follows: 



4^) = 2(a2)-^; af) = 2(«2) + ^, 
And the probability of success: 



(32) 



P?' = M [{af)' + (af + {if)' + (if )^) , (33) 

For the sake of simplicity, we can trace the effect of each iteration on the 
amplitudes of the system instead of writing the state of the system explicitly; 
for example, the amplitudes of the system after third iteration will be as 
follows: 

1- The mean of the amplitudes to be used in the diffusion operator: 



15 



(as) 



271+3 




-2- f ] 

V2 V 


^ TV y 



fe(2) 
V2 



(34) 



2- The new amplitudes: 



"o 


= 2 (as) - 


"o . 

72 ' 




(3) 
^2 


= 2 (as) - 


V2 ' 


af 


^(3) 


= 2 (ag) - 


"o . 


hf 


.(3) 
^2 


= 2 (ag) - 




^(3) 

"s 



(2) 

= 2(a3) + ^, 

(2) 

= 2(«3) + ^ 
2 (aa) 



j,{2) 

2(as) + ^. 



(35) 



3- And the probabihty of success: 

Pf=M((af)V(6f)');. = 0,1,2,3. (36) 

In general, after q iterations the recurrence relations representing the itera- 
tion can be written as follows: 

The initial conditions: o^q^ = b''^^ = 

1- The mean to be used in the diffusion operator: 

K) =^fi-^l;^>i- (37) 



2- The new amplitudes of the system: 



,(1) - O U \ I iq) _ o \ "o-.2g-2-l 



2(ai) + -^, a^!i2.-i-i = 2 (a,) T i ? > 2- (38) 



= 2 (a,) - ^, 6( = 2 («,) T g > 2. (39) 
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3- The probability of success for g > 1: 

Pi'^) = M [{a^^)' + (feS^))') ■t = 0,1, 2, ...,2^-1 - 1. (40) 

Using mathematical induction, we can prove that the probability of suc- 
cess after q iterations shown in Eqn. ()40p can take this form (Appendix B): 

Pi^^=(^-lj(l-^j +1, (41) 

4.2.1 Performance of Iterating the Algorithm 

i- The case where multiple instances of a match exist within the search 
space: Consider the following cases using Eqn. (j¥T|) : 

1- The case where M = N/2: the algorithm can find a solution 
with certainty after arbitrary number of iterations (one iteration 
is enough). 

2- The case where M > N/2: the probability of success is; for in- 
stances, at least 92.6% after the first iteration, 95.9% after second 
iteration and 97.2% after third iteration. 

3- For iterating the algorithm once {q = 1) and to get probability at 
least one-half, so, M must satisfy the condition M > N/8. 

ii- The case where few instances of a match exist within the search space: 

First, we need to represent the number of iterations q in terms of the 
ratio M/N. From Eqn. ()4Hl and using Taylor's expansion we get: 

pil) _ M 

For the cases where g > 1, the following conditions must be satisfied: 

n > 4 and 1 < M < N/8. (43) 
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It means that first iteration will cover approximately 87.5% of the problem 
with probability at least one-half; two iterations will cover approximately 92% 
and three iterations will cover 94%. It is easy to prove that the rate of the 
increase of the coverage range will decrease as number of iterations increases 
as shown in Fig.©. We can also see from Eqn. p^ and Eqn. pHj) that the 
algorithm needs O (N/M) iterations for > 4 and 1 < M < A^/8, which is 
similar to classical algorithms behaviour. It leads to a conclusion that first 
few iterations of the algorithm will do the best performance and there will 
be no big gain from continuing to iterate the algorithm. 




M/M 

Figure 7: A plot of the probability of success of the iterative version of the 
proposed algorithm where q = 1,2,. .. ,6. 



5 A Hybrid Quantum Search Engine 

We have devised a quantum search algorithm, which performs very well in 
case of multiple instances of the solution within the search space and a classi- 
cal behaviour in case of few instances of the solution. On contrary, Grover's 
algorithm needs O (^y^N/M^ to solve the problem but it's performance de- 
creases for M > N/2 [IHj. 
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This leads up to propose a hybrid quantum search engine, which combines 
both algorithms and can be integrated as follows: 

i- If the number of solutions M is known in advance: 

1- If 1 < M < N/8: Use Grover's algorithm with O (^^/N/Jiy 

2- If < M < N: Use the proposed algorithm with 0(1) . 

ii- If the number of solutions M is unknown: 

Iterate the proposed algorithm few times; say three iterations, which 
results in a chance of approximately 94% to find a solution. If it 
fails, we apply Grover's algorithm so we still have the same complexity 

O i^^/N0^. 




K A I 



Figure 8: The probability of success after five iterations from Grover's algo- 
rithm's vs. the proposed algorithm. 

We can see from Fig.® that Grover's algorithm is much faster in the 
case of few instances of the solution (ratio M /N is small) and the proposed 
algorithm is more stable and reliable in case of multiple instances of the 
solution. 
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6 Conclusion 



In this paper, we proposed a quantum search algorithm, which performs very 
fast in the case of muhiple instances of the solution within the search space 
(almost constant run-time) but the performance turns out to be classical for 
few instances of the solution. 

On the other hand we have Grover's algorithm, which performs very well 
in case of few instances of the solution and the performance decrease as 
number of solutions increase within the search space. 

This gave us the chance to propose a hybrid quantum search engine with 
general performance better that any pure classical or quantum search algo- 
rithm and still has O ^V^j for the hardest case and approximately 0(1) for 
M > N/8. 
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Appendix A 

To proof the identity shown in Eqn. (j26|) . we first need the next Lemma. 

Lemma 6.1 Let a and (3 angles; < < 7r/2, such that : 

sin^ [a) = cos^ . (44) 
Then, if k is any odd positive integer, we have: 

sin^ (ka) = cos^ (kp) . (45) 
Proof Since, < a, (3 < 7r/2, then we have, 

sin (a) = cos (/3) . (46) 

Also, we can write the following: cos (/5) = sin {ti/2 — j3) = sin (a). 

So, a = 71/2 — P . 

Therefore; 

cos (kp) = cos (f — «) ) 
= cos if - ka) 

= cos (^) cos (/ca) + sin (^) sin (to) ^ ' 

= ± sin (kp) . 

Then, cos^ {kp) = sin^ (A;/5) . 

Theorem 6.2 For any odd positive integer k and an angle O^i; < O^j < |, 
defined as follows: 

sin2(M = f- (48) 

Then, 

N 

^^CMsin2(A;M = 2"-^ (49) 



M=l 



Proof Consider 



sin" {Om) + sin^ (^^-m) = ^ + ^^^^ = 1- (50) 
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So, 

By Lemma, 
Or, 

Now consider, 



sin^ {9m) = 1 - sin^ (Oj^^m] 
= cos^ (On^m) ■ 



sin^ (kOM) = cos^ {k9 



N-M ■ 



sin^ (/c6'm) + sin^ {kON-jvi) = 1- 



N N N 

2 E ''CMsm^keM)= E ^CMsin2(A;^,,)+ ^ ^Cjv-m sin^ (A;^jv/ 

M=l M=l Af=l 

TV AT 

= E ^CMsin2(A;^,,)+ ^ ^Cm sin^ 

M=l M=Q 

N 

= E ^C^M (sin2 (A;^jvf ) + sin^ (A:^^„m)) 

M=l 

TV 

= E ^^M 
M=l 

= 2". 
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Appendix B 



To prove that the probabihty of success after q iterations is as shown in 
Eqn. (PT|) . we need first to prove the following relation: 

Let h^Q^ = = and given by the definition of the diffusion operator 
for g > 1 that, 

bl^^''^ ( M 



And, 



= 2 - ^ (56) 



Then, 



V2 



&r = 7^(l-2^) (57) 



Proof (By Mathematical Induction) 
Step 1: For g = 1, it follows directly from Eqn. (jlfi|l as follows: 

1.(0) 

= ^ (1-2^) 
Step 2: Assume the relation is true for q = t: 



(^) 

Step 3: Prove for g = t + 1: 
By definition, 



(58) 



C = 7^ 1-2^ (59) 



(60) 
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And, 



2!£ h _ ^^^ _ !£1 (61) 



V2 n) 

Substitute by the assumption, it directly gives the term for g = t + 1, 
e^^iri-2^V" (62) 



N 



Now, to prove that the probabihty of success of the proposed algorithm 
after q iterations can take this form: 



Given by definition that. 



Pi'^) = M {{a^^)" + ; ^ = 0, 1, 2, 2'^-' - 1. (64) 



Proof (By Mathematical Induction) 

Step 1: For g = 1, it is straight forward from Eqn.(24). 
Step 2: Assume the relation is true for q = t, 



(65) 



(#-l)(l-^) +1 



Step 3: Proof for q = t + 1, 
By definition. 
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Using Eqn. ljKTj) . we have, 

- ( h - (^ - 2A'n2* 

— 1 (^ _ M.\'^ (^ _ 2j\/^2t 
Substitute in Eqn.(|M|). we get, 

Pi-^) = M2-3_^ (1 _ Mf (1 _ + (f - 1) (1 - ^) 

= (f-i) (i-¥)"(^(f-i) + i) + i 
= (f-i)(i-¥)''''^ + i 
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